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Abstract 

Several recent results are reported from work aiming to improve the quantitative precision of relativistic viscous 
fluid dynamics for relativistic heavy-ion collisions. The dense matter created in such collisions expands in a highly 
anisotropic manner. Due to viscous effects this also renders the local momentum distribution anisotropic. Optimized 
hydrodynamic approaches account for these anisotropies already at leading order in a gradient expansion. Recently 
discovered exact solutions of the relativistic Boltzmann equation in anisotropically expanding systems provide a 
powerful testbed for such improved hydrodynamic approximations. We present the latest status of our quest for a 
formulation of relativistic viscous fluid dynamics that is optimized for applications to relativistic heavy-ion collisions. 
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1. Motivation 

Relativistic viscous hydrodynamics has become the 
workhorse of dynamical modeling of ultra-relativistic 
heavy-ion collisions. It is an effective macroscopic de¬ 
scription based on coarse-graining (via a gradient ex¬ 
pansion) of the underlying microscopic dynamics. Its 
systematic construction is still a matter of debate, com¬ 
plicated by the existence of a complex hierarchy of 
micro- and macroscopic time scales that are not well 
separated in relativistic heavy-ion collisions. Exact 
solutions of the highly nonlinear microscopic dynam¬ 
ics can serve as a testbed for macroscopic hydrody¬ 
namic approximation schemes. Such solutions have 
been found for the Boltzmann equation in the Re¬ 
laxation Time Approximation (RTA), which describes 
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weakly interacting systems, under the assumption of 
highly symmetric flow patterns and density distributions 
(Bjorken El and Gubser flows) of the system |[3H2|. 
Here we use both to test different hydrodynamic expan¬ 
sion schemes. 

2. Kinetic theory vs. hydrodynamics 

Hydrodynamics is an effective theory whose form is 
independent of the microscopic interaction strength. Its 
equations can thus be derived from kinetic theory in a 
window of weak coupling and small pressure gradients 
where both approaches are simultaneously valid. [We 
will here use the RTA Boltzmann equation as our start¬ 
ing point.] Only the values of the transport coefficients 
and the equation of state depend on the microscopic 
coupling strength; for the strongly coupled quark-gluon 
plasma created in heavy-ion collisions, they must be ob¬ 
tained with non-perturbative methods. 

The conserved macroscopic currents = (pf*) (par¬ 
ticle current) and 7’^*' = (energy-momentum 
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tensor) are obtained by taking momentum moments 
{h{p)) = f ^h{p)f{x, p) of the distribution func¬ 
tion f(x, p) (g is the degeneracy factor.) Hydrodynamic 
equations are obtained by splitting the distribution func¬ 
tion into a leading-order contribution fo, parametrized 
through macroscopic observables as iHini 


fo(x,p) = fo 


sJPp'^''(x)Pv-fi(xy 


(1) 


and a smaller first-order correction Sf (|b///ol « 1): 


fix, p) = foix, p) + 6fix, p). (2) 


In Eq. Q, ’EP^x) = u^^{x)u''ix) - <^{xW\x) + 
where the hydrodynamic flow field u^(x) defines the 
local fluid rest frame (LRF) and A'^'' = g^^''—u^u'' is 
the spatial projector in the LRF. <l)(jr) and the tensor 
fix) partially account for bulk viscous effects and shear- 
viscous deviations from local momentum isotropy in 
anisotropically expanding systems. T{x), p{x) are the 
effective temperature and chemical potential in the LRF. 

M^(x), f{x), and fi{x) are fixed by the Landau match¬ 
ing conditions, requiring = 1: 

Tfu'' = &uf {u-p)sf = {iu-p)\f = 0. (3) 


Here the eigenvalue S(T,p-, f, d)) is the LRF energy den¬ 
sity. The true local temperature 7’(7',/i; O) and local 
chemical potential are introduced by de¬ 

manding S(f,p;f,^) - SeqiT,p) and N{f,fl;f,<i) = 
{u-p)fy = Hoif, <l))A/eq(r,yu) where £eq, A4q are the ther¬ 
mal equilibrium energy and particle densities and 'Rq is 
a factor that depends on the viscous deformations f and 
O of the local momentum distribution ElIIollIIl. 

Writing = T^''+6Tf^'' = T^''+W'', / ^ = 

the conservation laws 


d^Ti^fx) = 0, 


d^fix) = 


Njx) - 7Veq(x) 

Trel(2c) 


(4) 


are sufficient to determine m^(x), T(x), pix), but not the 
dissipative corrections <1), H'^'', and V^', their evolu¬ 
tion is controlled by microscopic physics. Different hy¬ 
drodynamic approaches can be distinguished by the as¬ 
sumptions they make about the dissipative corrections 
and/or the approximations they use to derive their dy¬ 
namics from the underlying Boltzmann equation: 

1. Ideal hydrodynamics assumes local momentum 
isotropy of f, setting fo to be isotropic = 0) and 
all dissipative currents to zero: <h = H^*' -V^-0. 

2. Navier-Stokes (NS) theory maintains local momen¬ 
tum isotropy at leading order (i.e. in fo), sets <1) = 0, and 
postulates instantaneous constituent relations for H^*' 


and by introducing viscosity and heat conduction as 
transport coefficients that relate these flows to their driv¬ 
ing forces, but ignores the microscopic relaxation time 
that is needed for these flows to adjust to their Navier- 
Stokes values. This leads to acausal signal propagation. 

3. Israel-Stewart (IS) theory [TT] improves on NS the¬ 
ory by evolving H'^’' and y^ dynamically, with evo¬ 
lution equations derived from moments of the Boltz¬ 
mann equation, keeping only terms linear in the Knud- 
sen number Kn = Tnifp/Tmacro- 

4. Denicol-Niemi-Molnar-Rischke (DNMR) theory 

EH improves IS theory by keeping nonlinear terms up 
to order Kn^, Kn ■ Re ' when evolving H'^'" and y^. 
(Terms of second order in the inverse Reynolds number 
Re * vanish in the RTA used here but would otherwise 
appear, too.) 

5. Anisotropic hydrodynamics (aHydro) lfT3l allows 

allows for a leading-order local momentum anisotropy 
if^'’, <1) 0), evolved according to equations obtained 

from low-order moments of the Boltzmann equation, 
but ignores residual dissipative flows: YW'’ - -0. 

6. Viscous anisotropic hydrodynamics (vaHydro) lfT4l 
improves on aHydro by additionally evolving (using 
IS or DNMR theory) the residual dissipative flows 
nt'’', generated by the deviation 6f around the lo¬ 
cally anisotropic leading-order distribution function /q. 


3. Exact solutions of the Boltzmann equation 

3.1. Systems undergoing Bjorken flow 

For highly symmetric flow profiles the Boltzmann 
equation can be solved exactly in RTA. Bjorken flow 
El describes the dynamics of a longitudinally boost in¬ 
variant, transversally homogeneous system. A system 
with these symmetries is most conveniently discussed 
in Milne coordinates (T,r,tp,ri) where t = (f^-z^)'^^ 
and T] - i ln[(f-z)/(f-i-z)]. Bjorken flow is static in 
these coordinates, u^^ - (1,0,0,0). In Cartesian co¬ 
ordinates this implies a longitudinal flow velocity pro¬ 
file Vj = z/f El- The metric in Milne coordinates is 
ds^ = dr^-dr^ - r^dtp^ - T^drf. The Bjorken sym¬ 
metry restricts the possible dependence of the distri¬ 
bution function fix,p) to fix,p) - f{T\ Pi_,w) where 
w = tp^ - zE - Tmj_ sinh(y-; 7 ) is boost invariant. The 
RTA Boltzmann equation then simplifies to an ordinary 
differential equation 


< 5 t/(t;p±,w) = - 


fiT;pj_,w)-feqiT;pj_,w) 

Trel(T) 


with the solution ||3l|4l 


(5) 


fir, Px,w) 


= Dir, To)foipj_,w) (6) 

r dr' 

+ ( -;r77T>(T,T )/eq(T ;Pi,H’) 

Jto Ael(T') 
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where D(t 2 , ti) = exp(- J/) ^). 


3.2. Systems undergoing Gubser flow 

For systems with longitudinal boost invariance and 
azimuthally symmetric radial dependence, Gubser ^ 
found a flow pattern describing simultaneous longitudi¬ 
nal and transverse radial expansion by starting from a 
static system, - (1,0,0,0), in de Sitter coordinates 
(p, 0, (f), jf), with 


P(t, r) - — sinh 
0{t, r) - tan“* 


i/i-^V + ^V 


2qT 


2qr 


1 -H q^r"^ - q^r^ 


(7) 


In Cartesian coordinates this leads to the flow profile 
Vz - zfl, Vr - where q is an arbitrary scale 

parameter. The de Sitter metric is ds^ = ds^ - 
dp^- cosh^p {dO^+s,\r^ 0 d(p^)-dq^. A distribution func¬ 
tion that preserves the symmetries of this metric can 
depend on {x,p) only as f(x,p) - f{p',p^,pfl where 

Fn “ Fe + ^nd = w. Due to conformal symmetry 
of the metric and flow profile ||2l the particles described 
by / must be massless. With T(t, r) = T(p(t, r))/T the 
RTA Boltzmann equation simplifies to the ordinary dif¬ 
ferential equation 


dp 


fiP'^Pn^Pi) = - 


tip) 

c 


f{p'^Pn^p/)- U 


Tip)! 


where for the conformal case Trei(jc) = clT(x) and c/5 - 
77 /s is the specific shear viscosity. The solution is El 


fip;p^,w) = D(p,po)fo{p'^,w) ( 8 ) 

-I-- r dp't(p')D(p,p')fl^(p'-,pj^,w). 
c Jp„ 


This integral equation in p can be solved separately for 
each pair of momenta {p^,w). It can be worked out 
for any “initial” condition foipQ, w) = f(po', p^, w). We 
here use equilibrium initial conditions, /o = flq. 


3.3. Hydrodynamics of Gubser flow 

By taking hydrodynamic moments, the exact / yields 
the exact evolution of all components of Here, due 
to the high degree of symmetry, H^’' has only one in¬ 
dependent component, the longitudinal shear stress 
This exact evolution, which reflects the exact solution 
of the microscopic dynamics, can be compared to so¬ 
lutions, in Milne or de Sitter coordinates, of the vari¬ 
ous hydrodynamic approximation discussed at the end 
of Sec.[^ using identical initial conditions. For Bjorken 


flow, the equations corresponding to the hydrodynamic 
approximations 1.-6. are given and their solutions com¬ 
pared to the exact result in ||4l|5][T4l. For Gubser flow 
one finds the following El HU US: 

1. Ideal hydrodynamics: fideaiip) - fo/cosh^^^(p). 

2. NS theory solves instead the differential equation 
(1 /t)(df /dp) + |tanhp = l^p(p) tanhp for viscous 
T-evolution ||2l [15], with s n/j/its) and = 
5/7 tanhp where q/s- jTfrei (all variables with hats are 
made unitless by multiplying with appropriate powers 
of t). 

3. In IS theory the instantaneous constituent equation 
for in NS theory is replaced by a dynamical evolu¬ 
tion equation E): 

dn/^/dp + 5 {nj) tanhp + ^^/frei = ^ tanhp. 

4. DNMR theory adds to this evolution equation a term 
of order Kn-Re 

dn/^/dp+^{n/) tanhpH-Tr^/frei = tanhp-i-^7r[J tanhp. 

5. Anisotropic hydrodynamics (aHydro) solves a dif¬ 
ferent set of evolution equations for the temperature and 
pressure anisotropy which can be found in lIThl . 

6. For viscous anisotropic hydrodynamics (vaHydro) 
the form of the hydrodynamic equations corresponding 
to Gubser symmetry have not yet been worked out. 

4. Results 

For Bjorken flow the comparison of the exact evo¬ 
lution with the approximations 1.-6. was reported in 
mill] [Ml and will not be repeated here due to space re¬ 
strictions. It was found that NS theory reproduced the 
exact result better than ideal fluid dynamics, second- 
order viscous hydrodynamics (IS or DNMR theory) 
works better than first-order NS theory but failed badly 
for systems with very large viscosity (it leads to large 
entropy production even in the free-streaming, infi¬ 
nite viscosity limit), that anisotropic hydrodynamics re¬ 
paired this deficiency and provided qualitatively correct 
descriptions of the exact result for both small and large 
viscosities, but that all these approaches were outper¬ 
formed by vaHydro which, for Bjorken flow, reproduces 
the exact result almost perfectly over the entire range of 
viscosities and initial conditions m. 

Here we discuss one example from the analogous 
comparisons made in El fThll for systems expanding with 
the Gubser flow profile. The comparison, shown in 
Fig. is made in de Sitter coordinates - the translation 
to Minkowski coordinates is straightforward, and corre¬ 
sponding plots can be found in Ellll- Similar to what 
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47ir|/s = 1 



4jtr|/s = 3 




ent expansion rates. In fact, the transverse expansion of 
the Gubser solution is extremely rapid, with transverse 
velocity gradients that produce large dissipative flows 
and therefore severely stress any hydrodynamic frame¬ 
work based on an expansion around local thermal equi¬ 
librium, probably more severely than the flow profiles 
occurring in real heavy-ion collisions. That anisotropic 
hydrodynamics can capture these large dissipative flows 
and evolve them with high physical accuracy, as illus¬ 
trated by the excellent agreement with the exact Gub¬ 
ser solution, suggests that it is a robust framework that 
yields reliable results even under extreme conditions, 
such as those existing in the expanding quark-gluon 
plasma fireballs created in small collision systems such 
as p-i-A or p-Hp. 


Figure 1: de Sitter time (p) evolution of the effective temperature t 
(top row) and the scaled shear stress component (denoted by 
in the figure) (bottom row), for the exact solution of the Boltzmann 
equation (sold black line), Israel-Stewart (short-dashed blue line) and 
DNMR (dash-dotted green line) second-order viscous hydrodynam¬ 
ics, and anisotropic hydrodynamics (long-dashed red line). The two 
columns correspond to different specific shear viscosity values as in¬ 
dicated. Locally isotropic thermal equilibrium initial conditions with 
f = 0.002 were imposed at po = -10. (Taken from Ref. 1161 .) 

was reported above for Bjorken flow, anisotropic hydro¬ 
dynamics reproduces the exact microscopic Boltzmann 
dynamics much more accurately than standard second- 
order viscous hydrodynamics, especially for systems 
with large viscosity and pressure anisotropy (see my 
This happens even before the effects from residual dis¬ 
sipative flows, resulting from the deviation of the distri¬ 
bution function from the leading-order ansatz /o made 
in aHydro, are included. vaHydro includes the latter 
but has not yet been worked out for Gubser flow. As 
in the Bjorken case, it is expected to further improve 
agreement with the exact solution. 

A recent analysis of the exact solution of the Boltz¬ 
mann equation with Gubser symmetry showed that, for 
thermal equilibrium initial conditions at pQ, the solution 
becomes unphysical for large negative p-pu, as a result 
of the distribution function turning negative for large 
transverse and small longitudinal momenta Q. With 
such initial conditions it should therefore only be trusted 
for positive de Sitter times, p-po > 0 (as in Fig.[^. 

As seen in Fig.[^ where the exact solution is physical, 
it can serve as a valuable benchmark for assessing the 
accuracy and efficiency of various hydrodynamic ap¬ 
proximation schemes. The power of the Gubser solu¬ 
tion lies in the fact that it embodies a key feature of rel¬ 
ativistic heavy-ion collisions, namely simultaneous lon¬ 
gitudinal and transverse expansion with sharply differ¬ 


5. Conclusions 

Comparing various hydrodynamic approximation 

schemes with exact solutions of the Boltzmann equation 

we found the following performance hierarchy; 

vaHydro > aHydro > DNMR ~ IS > NS > ideal fluid 

We expect this to continue to hold for strongly coupled 

liquids which cannot be described by kinetic theory. 
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